The Multi-Dimensional Ensemble Empirical Mode Decomposition Method

A multi-dimensional ensemble empirical mode decomposition (MEEMD) for multidimensional data (such as images or solid with variable density) is proposed here. The decomposition is based on the applications of ensemble empirical mode decomposition (EEMD) to slices of data in each and every dimension involved. The final reconstruction of the corresponding intrinsic mode function (IMF) is based on a comparable minimal scale combination principle. For two-dimensional spatial data or images, f(x, y), we consider the data (or image) as a collection of one-dimensional series in both x-direction and y-direction. Each of the one-dimensional slices is decomposed through EEMD with the slice of the similar scale reconstructed in resulting two-dimensional pseudo-IMF-like components. This new two-dimensional data is further decomposed, but the data is considered as a collection of one-dimensional series in y-direction along locations in x-direction. In this way, we obtain a collection of two-dimensional components. These directly resulted components are further combined into a reduced set of final components based on a minimal-scale combination strategy. The approach for two-dimensional spatial data can be extended to multi-dimensional data. EEMD is applied in the first dimension, then in the second direction, and then in the third direction, etc., using the almost identical procedure as for the two-dimensional spatial data. A similar comparable minimal-scale combination strategy can be applied to combine all the directly resulted components into a small set of multi-dimensional final components. For multi-dimensional temporal-spatial data, EEMD is applied to time series of each spatial location to obtain IMF-like components of different time scales. All the ith IMF-like components of all the time series of all spatial locations are arranged to obtain ith temporal-spatial multi-dimensional IMF-like component. The same approach to the one used in temporal-spatial data decomposition is used to obtain the resulting